a^-dependence of J/ty suppression in pA collisions. 
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Abstract 

Coherence effects are important in the description of collisions with ex- 
tended objects as nuclei. At fixed target energies and small xp, the coherence 
length of the fluctuation containing the cc is small and the usual nuclear ab- 
sorption model is valid. However, at higher energies and/or xp the nucleus is 
seen as a whole by the fluctuation. In this case, the total, not the absorptive, 
cc — N cross section controls the suppression and also shadowing of gluons 
appears. We propose that the growth of the coherence length can explain 
the Xir-dependence of present experimental data. For this, we need a ratio of 
absorptive over total cc — N cross section of 0.2. 
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The J/^f suppression is one of the main signals for quark gluon plasma (QGP) 
formation The anomalous suppression observed by NA50 Collaboration is 
interpreted as produced by a deconfined state |§ though hadronic interpretations 
are possible |3J] . Whether NA50 data give or not a definite proof of the formation of 
a deconfined state is a topic of intense discussion. In this letter we study the non- 
anomalous suppression in the whole range of xf a subject that, though at first sight 
may seem solved, has some remaining open problems. This normal suppression is 
usually ascribed to multiple interaction of the produced pre-resonant cc state (color 
octet ||) with the surrounding nuclear matter. The picture is very simple and 
well known: a cc is created at some point z$ inside a nucleus in an octet state. In 
its travel through the nucleus, this state can interact at points z > z with other 
nucleons that will destroy it with a cross section a abs . The formula describing this 
nuclear absorption is, after integrating in z: 



where Ta{o) = J dzpA{b, z) is the profile function for nucleus A normalized to 1 and 
PaQ>, z) is the nuclear density, that we take from ref. 0. 

This formula describes well the observed pA data at midrapidities measured by 
NA38-NA51 |7| both for J/ty and ty' suppression, supporting the interpretation 
of the preresonant color octet state. A cross section of a abs = 6.5 ± 1.0 mb is 
obtained by the experimental collaboration. This cross section does not depend 
on the energy. For this analysis, data from OCu, OU and SU collisions were also 
included. The E866/NuSea Collaboration || has also measured J/\l/ suppression in 
pA collisions and the result does not fully agree with the one from NA38-NA51: using 
the parametrization a pA = A a a pp one obtains a ~ 0.95 for E866/NuSea || at y ~ 
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and a ~ 0.92 for CERN data. To reproduce E866 data at xf ~ a a abs ~ 3 mb is 
needed. E866/NuSea data give also the xf dependence of the nuclear suppression, 
the origin of which can not be completely attributed to modifications of the nuclear 
gluon distributions [f| (see bellow) and remains an open problem. The main goal of 
the present work is to describe this xf dependence and also to compare the different 
experiments. Understanding the rapidity pattern of the absorption is very important 
in order to have a real knowledge of the physics behind the suppression and also 
to extrapolate to RHIC and LHC energies, where this nuclear absorption would be 
present. 

The idea is simple and its theoretical formulation has been derived in a previous 
paper ||10|| . In the frame where the nucleus is at rest, the incoming proton fluctuates 
in a complicated system of quarks and gluons with coherence length l c . At small 
energies and xf, l c is small (of the order of the nucleon size) and only one nucleon 
in the nucleus takes part in the hard interaction that produces the cc. This implies 
<7pA ~ A. This behavior is modified by the collisions of the produced cc with the 
other nucleons in the nucleus. We have, in this way, the usual description of nuclear 
absorption given by eq. ([!]). However, at large energies and/or xp, l c gets eventually 
larger than the nuclear size and the nucleus is seen as a whole by the fluctuation. 
As a consequence the time ordering is lost and (|l|) is no longer valid. To describe 
this regime, we have introduced two types of collisions with the nucleons in the 
target, the ones of the light partons (mainly gluons) - with total cross section a - 
and the ones of the heavy system (cc) - with total cross section a - in an eikonal 
approach. The first ones give rise to modifications in the nuclear gluon distribution 
and the second to a suppression of charmonia states (the generalization of nuclear 
absorption). The result which replaces eq. (|I|) is [Tl| 
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= oSScFlbixi, Q 2 )9a(x2, Q\ b)e~^ b \ (2) 
with xp = Xi — x 2 , XiX 2 s = m 2 j^, and 

g A (x 2 , Q\ b) = 2 J duj [l - e -§<^)^Mb)j _ ( 3 ) 

uj represents kinematical variables of the gluon-nucleon interactions to be integrated. 
This term gives the modification of nuclear structure functions with respect to nu- 
cleons and will not be discussed here - see for instance |TT| for a model. In (|j) we 



see that the multiple scatterings of the heavy and light systems factorize, so that 
we can separate both contributions in the ratios R p a of pA to pp cross sections. 

D jjshadow rycc ( a \ 

tip A - ti v A K pA- I 4 J 

Neglecting shadowing corrections to gluons, the change from low to asymptotic 
energies consists in the substitution: 



1 



^abs 



1 _ exp A Ta(6) )] _ A Ta{p) exp (±j A TA{ p)\ } (5) 



2 

The most important point is the change of the absorptive cross section by the 
total one. When a = a abs the first correction term in the expansion in a is the 
same for both expressions. As these cross sections are not very large in practice, the 
numerical values turn out to be very similar. In this case, this justifies the use of 
formula ([I]) for high energy though strictly speaking it is only valid at small energies. 
In our case, however, we will suppose that a ^ a abs . That is, the preresonant cc 
state has a non negligible probability of not being destroyed. This increases the 
absorption for large values of xp where coherence is reached and a has to be used. 
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The two regimes are particular solutions of a more general equation which takes 
into account the coherence effects for any l c . In this case the factorization given in 
(U) is no longer valid. Neglecting shadowing to structure functions, 



^ = oj&Fgpfa, Q 2 )g P (x 2 , Q 2 ) £ C\ ± J <Pvr®{b)a® , (6) 



71=1 j=l ' 

where 

_ j_l — j 2 

rW) = f-otaNn-j 4- f-i - f--\ (_„abs\n-j+l 



a W = j j (-a**)*-' + (j - 1) (_ a ^)n-, + i . (7) 

The powers of nuclear profile functions in the expansion of @ have to be changed 
to 

r+oo p+oo r+oo n 

T n (b) — > Tj>\b)=n\l dzj dz 2 ... dz n cos(A(z 1 -z J ))Y[p A (b,z l ). 

J-oo Jzi Jzn-l j = i 

(8) 

The physical interpretation of the above equations is very clear, the first j — 1 
collisions are coherent (one of them corresponding to the light amplitude a), the 
last n — j are not and the j-th has the two possibilities. The A factor that controls 
these finite energy effects can be approximated by 

A.I=^ , (9) 

l c S Xi 

where m p is the mass of the proton and M 2 is the effective mass of the fluctuation. 
This effective mass could be measured in diffractive events containing a J/ty. In 
order to obtain the observed scaling of the suppression in xp, we take M 2 /s ~ 
constant. This means that increasing the energy allows the fluctuation to have bigger 
masses. This is the kind of behavior of triple pomeron contribution to diffraction 
Yl\ , however it is not clear that these are the relevant diagrams in the present case. 
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We take this scaling as a phenomenological ansatz, having in mind that this is most 
probably a finite energy effect that would disappear at asymptotic energies^. In 
practice, what we have done is taking M 2 = for the smallest experimental 

energy (-^/i ~ 20 GeV) and use the proportionality with s for the others. 

In our calculations we have also introduced the nuclear corrections to parton 
distributions (gluons) as given by the EKS98 parametrization ||. We could also 
use eq. (|]), however this would need some theoretical inputs from gluon-nucleon 
cross section that complicates the computation. Moreover, gluon antishadowing can 
not be reproduced by this formula. To take into account coherence effects we have 
defined 

(10) 

RpA S being the shadowing corrections to cc production computed with EKS98 
parametrizations, and A^ r j° b , A a e sy p and A^Jl the effective A (i.e. o v aIo vv ) com- 
puted with eq. ([I]), (fj) and respectively, all without shadowing. In this way, 

R shadow = 1 when ^ ig yery gmall ( A fin = A prob^ ^ R shadow = R EKS for 1&rge 

coherence lengths {A^ = A^" 1 ). However, as the data from E866/NuSea uses 
Be as a reference, the shadowing corrections are not large and, in any case, much 
smaller than the observed effect (see Fig. 1). The final result for absorption is given 
by eq. (§ with R* = A f e J}/A. 

With all these ingredients we have fitted the data from E866/NuSea Collabora- 
tion. The free parameters are a abs and a. To have a good description of E866/NuSea 

data we obtained <j abs = 3 mb and <?= 15 mb. The comparison with experimental 
■^1 thank A. Capella and A. Kaidalov for discussions on this point. 

6 



-pshadow 



a fin a asym 
^eff A eff 



A 



prob 
eff 



a asym 
A eff 



tdEKS 
n pA 



data can be seen in Fig. Also shown in this figure is the comparison taking 
into account only shadowing corrections to gluons given by EKS98 parametrization 
and nuclear suppression given by eq. ([!]) with cr abs =4.5 mb. Notice that if shadow- 
ing (antishadowing at xp ~ 0) is not included, o ahs ~ 3 mb is needed in order to 
reproduce the data at x^<0.2. 

In Fig. 2 we present the comparison with NA38-NA51 data. These data are 
measured in the interval 3 < yi a b < 4. We have used the values < xf >= 0.03 for 
Ei a b =450 GeV data and < xp >= 0.16 for E^ =200 GeV data. The description 
of the data is good in spite of the apparent discrepancy in the parameter a among 
the two sets discused above. Notice that the data at 200 GeV is more suppressed 
due to the larger xp. 

Finally, in Fig. 3 comparison is made with NA3 data |13| at Ei^ =200 GeV. The 
description is again not bad, though some discrepancy appears at large values of xp. 
This is also observed in other analysis []T4J and could be a signal of energy loss \ TE \, 
however, the evidence is too weak and more experimental data would be needed 
in this regionQ. Besides, including shadowing corrections (in fact antishadowing in 
this region) makes our results to increase, breaking the scaling on if mentioned 
above. This corrections are very uncertain and this could also be the origin of the 
discrepancy. 



Let's compare our analysis with previous ones. In |L4j a description of data 
very similar to ours was obtained by taking into account two different cross sections 
for octet and singlet cc states, the drawback of this model was the smallness of 
the lifetime of the octet state. In ||17||, some kind of energy loss of gluons from 



§ Actually, only W/Be data have been used in the fit, the Fe/Be data came out as a result. 
^Notice that energy loss with nuclear matter becomes smaller with increasing energy ]Tr| , so 
the effect in E866/NuSea will be smaller. 
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the incoming proton is proposed. In different effects as shadowing, energy 
loss, comover absorption and intrinsic charm are taken into account separately. In 
rP~gf| a model similar to ours, that takes into account coherence effects, has been 
developed. They also include two types of interactions for light and heavy systems, 
obtaining shadowing and cc suppression, however, their xp dependence come from 
loosing of momentum of the cc pair as it experiences multiple scattering. Finally, 



in p0[ coherence effects as well as other effects as energy loss and time formation 
are taken into account. This is a more formal analysis, but where coherence give 
effects similar to ours. Their computations are only for the \ c state. One main 
difference is that they don't take into account an effective mass of the fluctuation 
depending on the energy. This makes their comparison with E866 data less good 
than ours (we will obtain a similar result by fixing M 2 = Mj,^). In a previous work 
[ f21|j their comparison with E866 data where better but they limited their analysis 
to —0.1 < xp < 0.25 in order to see the influence of the formation time. 

To conclude, a high energy pA collision can be seen, in the laboratory frame, as 
the multiple scattering of the system of quarks and gluons into which the incoming 
proton fluctuates. The coherence length of this fluctuation increases with the energy 
and xf and two regions can be distinguished. When l c -C Ra no shadowing to gluons 
is present and the suppression of cc pairs is given by the usual formula ([l|). When 
I c >Ra the whole nucleus takes part in the collision and typical coherent phenomena 
as shadowing appear. The main point is that this change of regime is accompanied 
by a change a abs — > a in the expressions. If these two cross sections are different, the 
suppression at large values of xp and/or larger energies is bigger. Let's also comment 
about the experimental data from NA38 don't see any difference between J/\I/ 
and suppression, whereas some difference seems to appear in E866/NuSea. This 
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difference would be very easily accounted for in our approach just by taking the 
absorptive cross sections different for singlet and octet cc states: as the total J/ty 
has a contribution of ~40% singlet coming from disintegrations of Xc a larger a abs 
for octet than for singlet would explain the difference. This would introduce a new 
parameter. Finally, the extrapolations to RHIC and LHC energies depend on two 
assumptions, the energy dependence of a and a abs (that we have taken as constants) 
and the value M 2 of the effective mass of the fluctuation. With this last assumption 
we are wondering if the xf scaling of the suppression will still be valid at high 
energy (in the case M 2 / s ~ const.) or if, on the contrary, a scaling in x 2 will 
appear, this last possibility is the most reasonable as the proportionality in s of M 2 
seems to be a finite energy effect. Let's give some estimations for RHIC and LHC 
at Xp ~ 0. Assuming no energy dependence for the cross sections a and a abs we 
obtain for pAu collisions at y/s =200 GeV a ratio over pp of 0.81 if the scaling in 
x F is maintained and 0.43 if it is not. This means ratios of 0.67 and 0.17 for AuAu 
collisions respectively. In the case of LHC, the only difference is shadowing, that 
will not affect the case M 2 /s ~ const, in the other case, we obtain 0.31 for pPb 
collisions and 0.1 for PbPb at yfe =5500 GeV. 
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Figure captions: 



Fig. 1. Comparison of our results on xp dependence of nuclear suppression (solid 
lines) with E866 data || for Fe/Be and W/Be. Effects of shadowing as given by 
EKS98 parametrization including absorption with a abs = 4.5 mb in eq. ([[]) are also 
shown (dashed lines). 

Fig. 2. B^/A measured by NA38 Collaboration at E lab =A50 GeV 2 (black 
boxes) and Ei ab =200 GeV 2 (white boxes) compared with our results at i^ a6 =450 
GeV 2 (solid line) and ^^=200 GeV 2 (dotted line). Also shown, the suppression 
given by eq. (U) with a abs =3 mb. 

Fig. 3. a parameter as a function of xp compared with data from NA3 [T3| at 
£ /ab =200 GeV 2 . 
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Figure 1 
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Figure 2 
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Figure 3 
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